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Abstract 

We discuss pair production and decays of stops, sbottoms, r-sneutrinos, 
and staus in e + e~ annihilation in the energy range yfs = 500 GeV to 2 TeV. 
Numerical predictions within the Minimal Supersymmetric Standard Model 
(MSSM) for cross sections and decay rates are presented. We study the stop 
discovery potential for yfs = 500 GeV and 10 fb -1 integrated luminosity for 
polarized e~ beams. Moreover, we give an estimate of the error of the soft- 
breaking stop and sbottom parameters that can be obtained by cross section 
measurements with polarized e~ beams. 
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1 Introduction 


If nature is supersymmetric at the weak interaction scale, the masses of the supersym¬ 
metric (SUSY) particles are expected to be lower than approximately 1 TeV (I. |2J. 
The weakly interacting SUSY particles are then within the reach of an e + e - linear 
collider with a center-of-mass energy between 500 GeV and 2 TeV. An e + e~ linear 
collider in this energy range will not only be a discovery machine for SUSY particles 
i HQ- but will also allow detailed measurements of the underlying SUSY parame¬ 
ters H {?;. ||. 9- The experimental search of SUSY particles which are relatively light 
will be particularly important. The lighter scalar top quark U, the SUSY partner of 
the top quark, and for tan/3>10 also the sbottom b\ or the stau f\ may even be 
the lightest visible SUSY particle m m 0 [13]. The reason is that the Yukawa 


interactions reduce the soft SUSY breaking masses of the left and right sfcrmions, 
f L and fn, of the 3 rd generation, compared to those of the 1 st and 2 nd generation 
|b§ |U|, and also induce a mixing which may make one mass-eigenstate rather light. 
The production cross sections and the decay rates, and thus the discovery reach of 
these sfermions show a distinct dependence on the /i,-/# mixing angles 0, |16|]. If the 
gluino is heavier than t\ and b\ the most important decay modes of these squarks are 
those into quarks and neutralinos or charginos. 

We present results for the production of stops, sbottoms, r-sneutrinos, and staus 
in e + e _ annihilation at energies between /s = 500 GeV and 2 TeV. We also dis¬ 
cuss in detail the decays of these particles. Furthermore, we present an example for 
signal selection and background rejection for stop production at a/s = 500 GeV and 
C = 10 fb -1 . In addition, we investigate the possibility of determining the masses 
and mixing angle of stops. If SUSY particles are experimentally discovered, the mea¬ 
surement of their properties will be the most important step further. Polarization of 
the e~ beam plays an important role. As we will show, from measurements of the 
production cross sections with polarized e~ beams we can determine the masses and 
mixing angles of the stops and sbottoms and in turn the underlying soft SUSY break¬ 
ing parameters with good precision. Knowing the latter we will be able to test the 
theoretical hypotheses about the SUSY breaking mechanism. 

We perform our calculations in the framework of the Minimal Supersymmetric 
Standard Model (MSSM) Jl], |J. This contains the Standard Model (SM) particles, 
sleptons £ ± , sneutrinos £/ squarks q, gluinos g, two pairs of charginos xt , i — 1, 2, four 
neutralinos x° k , k = 1,... ,4, and five Higgs particles h°, H°, A °, H ± [|T1 [2] [TtTI . The 
phenomenology of stops, sbottoms, staus, and r-sneutrinos, and their decay products 
is determined by the following parameters: 


- the soft-breaking parameters My, My, Mq, M y, My, A T , A t , At,, which determine 

the mass matrices of the stau, stop, and sbottom systems, and the mass of the 
r-sneutrino, 

- the (soft-breaking) SU{2) and U( 1) gaugino masses M and M', 
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the higgsino mass parameter /i, and tan f3 = 1 ’ 2 /vi (where V\ and v? are the vacuum 
expectation values of the neutral members of the two Higgs doublets). 


We assume the GUT relations M'/M = |tan 2 @jy ~ 0.5, and rrig/M = a s /oi 2 ~ 3, 
where trig is the gluino mass. Furthermore, we assume that the Xi is the lightest SUSY 
particle (LSP). 

The lower model independent mass bound for stops obtained at LEP1 is 45 GeV 
18fl . Stronger limits up to 65 GeV are reported from the data taking at LEP1.5 at 


130-161 GeV [|KJ]. The D0 experiment at the TEVATRON excludes the mass range 


40 GeV < mj< 100 GeV for the stop, if the mass difference rrif — m^o > 30 GeV [[20 


I 11 Section 2 we shortly review the basic facts about L-R mixing of stops, sbottoms, 
and staus. We also present numerical results for the production cross sections of 
stops, sbottoms, r-sneutrinos, and staus for unpolarized and polarized e - beams. 
In Section 3 we describe the decays of stops, sbottoms, r-sneutrinos, and staus and 
present numerical results for the most important branching ratios. We also list the 
signatures which are expected to be relevant at y/s = 500 GeV. In Section 4 we 
describe an event generator for t\ t\ production and decay. In Section 5 experimental 
sensitivities are determined based on Monte Carlo simulations. In Section 6 we give 
an estimate of the experimental error to be expected for stop and sbottom masses and 
the stop mixing angle, as well as for the soft SUSY breaking parameters. Section 7 
gives a summary. 


2 Cross Sections for Pair Production of Stops, 
Sbottoms, r-Sneutrinos, and Staus 


2.1 Left-Right Sfermion Mixing 


The SUSY partners of the SM fermions with left and right helicity are the left and 
right sfermions. I 11 the case of stop, sbottom, and stau the left and right states are in 
general mixed. I 11 the (//,, /^) basis the mass matrix is pll O 


M) = 


2 

m7 

JL 


afmf 


a f m f m } R 


( 1 ) 


with 


m 2 j L = M| + m| cos 2/3(Tj — e/ sin 2 6 W ) + m 2 , (2) 

rnj R = Mp, + e/m| cos 2/3 sin 2 9 W + mj , (3) 

a t = A t — /i cot /3, db = Ab — p tan /3, a T = A r — ptan/3, (4) 


where e/ and Tj are the charge and the third component of the weak isospin of the 
sfermion /, Mp = Mq for f L = t L , b L , Mp = M~ L for f L = f L , Mp, = Mp, Mp, Mp 
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for fn = in, bn, ?r, respectively, and mj is the mass of the corresponding fermion. 
From renormalization group equations |15] one expects that due to the Yukawa in¬ 
teractions the soft SUSY breaking masses Mq, M p, Mq, Mp, and Mq of the 3 rd 
generation sfermions are smaller than those of the 1 st and 2 nd generation. Evidently, 
tj r tn mixing can be important because of the large top quark mass. For sbottoms 
and staus L-R mixing can be important if tan/3 > 10. The mass eigenvalues for the 
sfermions / = t, b, t are 


m /i. 2 = 5 ( m h + =F 




1 }l 


1 )r 


~ m / i? ) 2 + 4m 2 a 2 ) 


(5) 


where t\, b\ and fi denote the lighter eigenstates. The mixing angles 0j are given by 

^,2 


cos 9j = 


if m? < m? . 

JR JL 


-a f m f 


(' m 2 jL — + a f m f 


1_ 

f L fl 

(6) 

/ 

sJ( m2 f L ~ m2 f j) 2 + a / m / 

1 

' 75 

if m? < rrif and cos#? 

JL JR 1 / 1 

< 75 


The mass of the r-sneutrino is 




ml T = Ml + \m 2 z cos 2 (3. 


(V 

srsions 

of 

eqs. (|D to (IP are 



M h 

= 

m 2 2 sin 2 6f + cos 2 9- t — m 2 z cos 2/3 — | sin 2 0 w ^j 

-m 2 t 

(8) 

Ml 

= 

ml sin 2 6i + m 2 2 cos 2 # f - — | m 2 z cos 2/3 sin 2 6\y — m 2 


(9) 

M l 

= 

rri~ 2 sin 2 0 b + rrij n cos 2 0 b + m 2 z cos 2(3 Q — | sin 2 0 w ^j 

-m 2 b 

(10) 

Ml 

= 

sin 2 0 b + m~ 2 cos 2 0 b + | m 2 z cos 2(3 sin 2 0 W — ml 


(11) 

Ml 

= 

m~ 2 sin 2 Of + m 2 fl cos 2 9 f + m 2 z cos 2(3 — sin 2 0 w ^j - 

- ml 

(12) 

Ml 

= 

m fl sin Of + rrif 2 cos 0 f + m z cos 2(3 sin 0 W — m T 


(13) 

Ml 

= 

m l T ~ 1 m i cos 2/3 


(14) 

a f m f 

= 

(m ^ — m 2 j 2 ) sin 0j cos 0j 


(15) 


The soft-breaking parameter Mq enters the equations of mi L and m bh . Therefore eqs. 
(§) and (|I0| ) imply the following condition (at tree-level): 

m 2 v cos 2/3 = sin 2 9~ t + cos 2 9~ t — mj- sin 2 6~ b — rnf^ cos 2 0 b — m 2 t + m 2 . (16) 

This condition shows that if tan (3 and five of the six measurable parameters 
rnp, mf 2 , cos 2 0i, m bi ,m b<2 , and cos 2 0 b are known, the sixth can be predicted. An anal¬ 
ogous condition holds for Mp in the slepton sector due to eqs. flT2|) and (|14|). Further¬ 
more, from eqs. © and (|15|) one can see that, in the convention used, a/ and cos 0j 
have opposite signs. 
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2.2 Cross Sections for e + e' 


fi fi 


The reaction e + e —> f t f t proceeds via 7 and Z exchange. The tree-level cross section 
at a center-of-mass energy y/s is given by |2T|, [ 2 / 


gtree _ 


na 2 Nc 
3 s 


P 3 


Q) + 


\vl + al)v\ 2 

if:, 4 f 4 S 

1 U 


Qf V e Vj 

2 Stf/Ciy 


■sis 


171 z) 


(s — m |) 2 + T|m| 


(17) 


where = 1 — c^y = sin 2 9w, v e = 2 sin 2 6w — \, and a e = — / IVc is a colour 
factor which is 3 for squarks and 1 for slcptons. The tree-level cross section has 
the typical /3 3 kinematic suppression where (3 — (1 — Amp/s) 1 ^ 2 is the velocity of 

the outgoing scalar particles. The Z coupling to f l f l is proportional Vj. with vj i = 
2 (Ij cos 2 9j — Qf sin 2 9w), vj 2 = 2 (I'j sin 2 Oj — Qj sin 2 9w), where I'j and Qf are the 
third component of the weak isospin and the charge of the fermion / (Q e = —1). This 
coupling vanishes at the mixing angles cos 2 9j = 0.31, 0.16, and 0.46 for f 1; 6 1; and 77, 
respectively. Note that the sign of cos Of cannot be determined from the cross section 
(O), as this depends only on cos 2 6j. 

The interference between the 7 and Z exchange contributions leads to a character¬ 
istic minimum of the cross sections for e + e~ —> which occurs at a specific value 

of the mixing angles 0j given by 

cos 2 9j\ min = sin2 Qw[ 1 + (1 - s/m|)F(sin 2 9 W )\ (18) 


where F(sin 2 #yv) = cos 2 6 w (L e + R e )/(L 2 e + R 2 e ) ~ —0.22, F(sin 2 6 W ) = cos 2 9 W /L e ^ 
—2.9, and F{sin 2 6w) = cos 2 9w/R e ~ 3.3, for unpolarized, left- and right-polarized 
e~ beams, respectively, with L e = —\ + sin 2 0w and R e = sin 2 9\y. For polarized 
e~ beams the dependence on the mixing angles is much more pronounced than for 
unpolarized beams. The corresponding minimum of the e + e _ —> /d,/d cross sections 
occurs at 1 — cos 2 9j\ min . 

I 11 the calculations of the cross sections we have also included SUSY-QCD cor¬ 
rections taking the formulae of ^3[ (see also pi], |24], ^5|) with a s (M z ) = 0 . 12 , 
m t = 175 GeV, and corrections due to initial state radiation H?7f| . 

Figure 1 a shows contour lines of the total e + e“ —> t\ti cross section in the 
mi — cos 9i plane for yfs = 500 GeV and unpolarized beams. For the calculation 
of the SUSY-QCD radiative corrections we have assumed mf 2 = mg = 300 GeV. Sig¬ 
nificantly above the threshold there is a clear dependence on cos Of. Figure lb shows 
the cos 9i dependence of the e + e“ —> t\i\ cross section for left- and right-polarized 
as well as for unpolarized e~ beams for y/s = 500 GeV and m .7 = 180 GeV. For 
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both left- and right-polarized e~ beams the cross sections depend strongly on the 
mixing angle. It is important to note that this dependence is opposite for left and 
right polarization. Therefore, experiments with polarized e~ beams will allow a more 
precise determination of the mass and the mixing angle Of. 

Figure 2 a shows contour lines of the total cross section of e + e“ —> f 2 ^2 in the 
rrif 2 — cos 2 Of plane at a/ s = 2 TeV. The cos Of dependence of this cross section at a/s = 
2 TeV for left- and right-polarized and unpolarized e~ beams is shown in Fig. 2 b for 
mf 2 = 700 GeV. Again, the cos Of dependence is much stronger for polarized than for 
unpolarized beams, however, the behavior is opposite to that of e + e _ —> t\t\. For the 
calculation of the SUSY-QCD radiative corrections we assumed = 300 GeV and 
rrig = 700 GeV. 


We shall discuss in Section 6 how experimental data on cross sections with left— 
and right-polarized beams would allow to determine masses and mixing angles, and 
then give information on soft SUSY breaking parameters. 

Figure 3 a shows the contour plot of the total cross section of e + e _ —> bibi in the 
rrif n — cos 0 h plane at /s = 500 GeV for unpolarized beams. For a polarized e~ beam 
the cos Of dependence of the cross sections is much stronger, as shown in Fig. 3 b for 
m~ bi = 180 GeV and a/s = 500 GeV. 


Contour lines for the cross section of e + e 


t\T i in the — cos Of plane 


at a/s = 500 GeV are shown in Fig. 4 a. Figure 4 b shows the cross section for 
polarized and unpolarized e~ beams as a function of cos Of for = 180 GeV and 
a/s = 500 GeV. For both beam polarizations these cross sections again exhibit a 
strong dependence on the mixing angle. 

Figure 5 shows the mo dependence of the cross section for e + e _ —> v T v T for 
unpolarized as well as left- and right-polarized e~ beams. Since the v r is not mixed, 
the polarization dependence is entirely due to the different Ze + e~ couplings. 

The a/s dependence of the e + e _ —> UU cross section is shown in Fig. 6 for = 
180 GeV and cos Of = 0.7. The effect of SUSY-QCD corrections for e + e~ —*• 
from gluon and gluino exchange as well as the initial state radiation correction as 
a function of /s is demonstrated in Fig. 7, for cos Of = 0.7, = 180 GeV, and 


m - 


t2 


= = 300 GeV. Note that at high energies the gluino exchange contribution 


has the opposite sign of the gluon exchange contribution, and the absolute values are 
increasing with /s. The effects are similar for e + e~ —> / b\. A detailed discussion of 
SUSY-QCD corrections is given in p3|. The effect due to initial state radiation turns 
out to be of the order of 10 


The sum of all corrections can well exceed 10%. 
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3 Decays of Stop, Sbottom, r-Sneutrino, and Stau 

The sfermions of the third generation can have the weak decays (i, j = 1, 2; k = 

1 _4 ) 


ti 

-»■ txl> 

bxt 

(19) 

bi 

- bxl 

t xj 

( 20 ) 

Ti 

-»■ T xt 

V rXt 

( 21 ) 

U T 

-»■ v x°k> 

TXj 

( 22 ) 


Owing to the Yukawa terms and the L-R mixing the decay patterns of stops, sbottoms, 
and staus will be different from those of the sfermions of the first two generations 
27], |28|. Stops and sbottoms can also have the strong decays 


tg, 


bg. 


(23) 


They are dominant if they are kinematically allowed. Otherwise, the lighter squark 
mass eigenstates decay mostly according to ([L9]) and (|20|) . Moreover, in case of strong 
L-R mixing the splitting between the two mass eigenstates may be so large that the 
following additional decay modes are present 0 : 


h ~ 

t x Z{h\ H°, A 0 ), 

h IU + ( H + ), 

(24) 

b 2 - 

- bi Z (h°, H°, A 0 ), 

ti W- (H-), 

(25) 

t 2 - 

-* h Z, 

v T W~. 

(26) 


If the t\ is the lightest charged SUSY particle and m^o +mb + mw < wy, < m^o + m t , 
the decay t\ —> bW + Xi can be dominant |29|, otherwise the higher-order decay 
ti —> cxi dominates p2| . In the case that rrif i < mg also t\ —> bu T f\ and in case 
that me, < mg also t\ —> br + i> T might be important. 

Figure 8 a and 8 b show the parameter domains in the M — g plane for the decays 
of t\ and b\, eqs. ([L9]) and ([ 20 ]), taking mg = 180 GeV, tan/3 = 2 , and mg = 180 GeV, 
tan/3 = 30. In region (a) only the decay U —> cxi is allowed, whereas in region (b) 
f i —> cxi and U —> b W + x® are possible. In the small stripe of region (c) also t\ —> 
is possible. In region (d) the decay t\ —> bxt l ias practically 100% branching ratio. 
This is further illustrated in Fig. 9 a where we plot the branching ratios for t\ decays 
as a function of M for /i = —500 GeV, mq = 180 GeV, tan/3 = 2, and cos 6 ^- = 0.7. 
The parameter domains for the f\ decays into neutralinos are almost identical to 
those of the corresponding bi decays, if the masses of f\ and b± are equal. 

Figure 9 b shows the branching ratios for bi decays as a function of M for rrij n = 
180 GeV, tan/3 = 30, cos 6 *^ = 0.7, and g = —500 GeV. Figure 10a and 10b show 
the branching ratios for b\ and f\ decays as a function of M for g = —130 GeV 
and m~ b] = rrif x = 180 GeV, tan/3 = 30, cos 9 b = cos9f = 0.7. Similarly, Figure 10c 
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shows the branching ratios for u T decays taking mo T = 180 GeV, tan (3 = 2 , and 
H = —130 GeV. The LSP decays of b { and f\ are about 60% and 40%, respectively. 
The branching ratio for the visible decays of the v T is between 40% and 90%. 

The decay patterns of the heavier sfermion mass-eigenstates can be quite com¬ 
plicated, because all the decay modes of eqs. to (|26|) can occur. We calculate 
the different decay widths with the formulae of Refs. )T3j. |l(| |29[. Figure 11 a and 


lib show the branching ratios for t 2 decays as a function of M for m p, = 700 GeV, 
taking tan f3 = 2 , p = —1000 GeV, rrq = 150 GeV, Mq = 607 GeV, Mp = 360 GeV, 
Mfj = 850 GeV, and A t = 500 GeV. The masses and mixing parameters of the 
other particles involved then are mp = 258 GeV, cos 6 ^ = —0.468, = 608 GeV, 


m 


l 2 = 851 GeV, cos 9~ b = —0.999, m^o = 98 GeV, mu o = 165 GeV, m H ± = 169 GeV, 
cos a = 0.535. We have included radiative corrections to the Higgs masses and the 
Higgs mixing angle according to [|3(]] . For M < 200 GeV the decay t 2 —■> tg dominates, 
whereas for M > 200 GeV the decay t 2 —■► Zt\ has the largest branching fraction. It is 
interesting to note that in this case the decays t 2 —■► A 0 t\ and t 2 —> H° t\ have larger 
branching ratios than the decay into bxt which is the most important decay mode 
of t\. 


Figure 12 a and 12 b show the branching ratios of the b 2 decays for m b2 = 700 GeV, 
tan /3 = 30, p = —1000 GeV, = 150 GeV, Mq = 637 GeV, Mp = 320 GeV, Mp, = 
450 GeV, and A t = 500 GeV. The masses and mixing parameters of the other particles 
involved then are m bi = 350 GeV, cos 6~ b = —0.469, mp = 324 GeV, m t - 2 = 678 GeV, 
cos Of = —0.273, m/p = 118 GeV, mu o = 146 GeV, m^± = 176 GeV, and cos a = 
0.998. In this example the transition b 2 —> H~ t\ is dominant for M>200 GeV, and 
b 2 Z b\ and b 2 —>■ W~ t\ have larger branching ratios than b 2 —»t xt or ^2 — 1 ► bxi 2 - 
The branching ratios for the decays of the sleptons t 2 and v T are shown in Figs. 13 
a, b, and Figs. 14 a, b, for rrif 2 = 700 GeV and m„ T = 687 GeV, respectively, taking 
tan f3 = 30, p = —1000 GeV, tua = 150 GeV, Mp = 690 GeV, Mp = 490 GeV, and 
A t = 500 GeV. The mass of the lighter stau then is m fl = 480 GeV, and the mixing 
angle is cos Of = —0.213. The Higgs boson masses are the same as in Fig. 12. In these 
examples, the t 2 and v T decays into charginos and neutralinos have large branching 
ratios, and the decays t 2 —> u Xi and v T —> rxt dominate. The decay f 2 —>Zfi 
has a branching ratio of 5% to 10%, and h T —> rxt has a branching ratio of 35% to 
55%. It is interesting to note that in this example the branching ratio of the invisible 
decays mode v T —> v T Xi is always less then about 18%. 

Quite generally, the decay widths of the U, bi, and f \ decays into neutralinos and 
charginos, eq. ( 0 ) to (0). depend also on the sign of cos 6j. It may therefore be 
possible to determine the sign of the mixing angle by studying these decays. 

Table 1 lists the most important signatures for H, %, /%, and % for y/s = 500 GeV. 
If the decays t] —> b xt or f\ — > v T xZ occur, the xt w iH most probably be discovered 
first and thus its mass and couplings will be known. The decay t\ —> b W + x i leads to 
the same final states as i\ —> bxt (provided xt Xi i s n °l allowed). The decay 

v T —^ v>Xi is invisible. Thus, one-sided events can occur where one v T decays invisibly 




and the other one decays visibly into one of the final states given in table 1. 


Channel 

Signatures 

ti-*b xt 

1 6-jet + 1 £ + + t> T , 1 6-jet + 2 jets + 

t\ -> cx i 

1 jet + i> T 

b\ -► bXi 

1 6-jet + j> T 

bi^b x° 2 

1 6-jet + + ]/) T , 1 6-jet + 2 jets + f> T 

Tl^T Xl 

T + i> T 

n -*• t X 2 

t + + ]/) T , t + 2 jets + t> T 

rf -»• Vt Xi 

t~ + $ T , 2 jets + j> T 

"r -> T~ Xt 

T~ + £ + + t> Tl T~ + 2 jets + j) T 

Vt -> V X2 

£ + £~ + j) T , 2 jets + j) T , 7 + 


Table 1: Expected signatures for t\, b\, z> r , and fi production for 
y/s = 500 GeV. Owing to pair production all combinations 
of the corresponding signatures may occur. 


4 Stop Event Generation 


In this section we describe the event generator for e + e~ —> t\ti with the stop decay 
modes t\ —> and t\ —> bxf. The chargino decays via xt W + x±, where W + 
can be either virtual or real. The event generator is based on the calculation of 
the 4-momenta distributions of the stop decay products Xi c X?c and xtb Xib- The 
large effects of QCD corrections are included in the cross section calculation. Stop 
production and decay have been defined as new processes in the PYTHIA program 

0 


package 

states. 


The event generation process includes the modelling of hadronic final 


In the first step of the event generation, initial state photons are emitted using 


the program package REMT |H] which takes into account the expected stop cross 
section from zero to the nominal center-of-mass energy. Beamstrahlung photons are 
generated using the beam parameters of the NLC 1992 design. The effective center-of- 
mass energy is calculated for the initial production of the 4-momenta of the final-state 
particles. These 4-momenta are then boosted to the lab-frame according to the mo¬ 
mentum of the emitted photons. For the hadronization process of the cc in the Xi c X?c 
and of the 66 in the xfb X\b decay mode, a color string with invariant mass of the 
quark-antiquark-system is defined. The possible gluon emission and hadronization 
are performed using the Lund model of string fragmentation with the PYTHIA pro¬ 
gram package j3l| . The Peterson et al. |52| fragmentation parameters for the c- and 
6-quarks are used: e c = 0.03 and = 0.0035. Finally, short-lived particles decay into 
their observable final state. Details of the event generator are given in |33] . 
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5 Simulation and Selection 


The investigated background reactions and their cross sections are shown in Fig. 15. 
They are simulated for £ = 10 fb -1 , and 1000 signal events are simulated in the Xi c X?c 
and xtbxib decay channels. The L3 detector at CERN including the upgrades for 
LEP2 served as an example for an e + e _ 500 GeV detector. Details of the parametric 

An important feature is the overall hadronic 


34 


detector simulation are given in 
energy resolution of about 7%. 

In both channels, the y°’s escape the detector and cause large missing energy. In 
the case of Xi c Xi^ the c-quarks form mostly two acoplanar jets. A mass combination 
of = 180 GeV and rn^o = 100 GeV is investigated in detail. For xtbxib on 
average the visible energy is larger. In this channel, the mass combination = 
180 GeV, m^+ = 150 GeV, and m^o = 60 GeV has been studied. Typically four jets 
are formed, two from the 6 -quarks, and two from the boosted IF’s. 

In the first step of the event selection, unbalanced hadronic events are selected 
using the following selection requirements: 


25 < hadronic clusters < 110, 0.2 < E vis /\fs < 0.7, 

Ef b / E yis < 0.5, Thrust < 0.95, | cos thrust I < 0.7 . 


Channel 

Xi c Xi c 

XibXi b 

qq 

WW 

eWi/ 

tt 

ZZ 

eeZ 

Total (in 1000) 

1 

1 

125 

70 

50 

7 

6 

60 

After preselection (in 1000) 

0.4 

0.7 

1.7 

2.2 

3.2 

1.3 

0.2 

0.3 


Table 2: Expected events per 10 fb 1 at -y/s = 500 GeV, and number of 
events after the preselection as defined in the text. 


A large part of the background of back-to-back events without missing energy is re¬ 
jected. Table |] shows the number of initially produced events per £ = 10 fb -1 at 
a/s = 500 GeV, and the number of events which pass this preselection. The require¬ 
ment of a large number of hadronic clusters removes e + e _ , /i + /i _ , and most of the 
t + t~ events. The minimum energy cut reduces most of the 77 events and ensures 
almost 100% trigger efficiency. The background from 77 events can, in addition, be 
strongly reduced by rejecting events where a scattered initial electron is detected at 
low angles. The upper energy cut reduces all standard background reactions. Beam 
gas events and events where much energy goes undetected along the beam axis are 
removed by rejection of events with very large parallel imbalance. The thrust cut 
removes remaining r + r~ events and reduces largely qq and Z°Z° background. The 
cos ^Thrust cut removes events where most probably much energy escapes undetected 
along the beam axis. 

The final Xi c XiC event selection is summarized in Table [|. The following cuts are 
applied: 
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• A hard upper energy cut reduces all standard background except eWv (Fig. 16). 

• Jets are clustered using the JADE algorithm. The y-cut value is optimized to 
obtain two jets for the signal. 

• Semileptonic decays of the top quark can induce missing energy. These events 
are partly removed by requiring no isolated electron or muon. 

• Events with large longitudinal energy imbalance are removed where probably 
much energy escapes undetected along the beam axis. 

• The invariant mass of the two jets is required to be larger than 120 GeV to 
remove almost entirely eWv events (Fig. 17). 

• The acoplanarity angle is defined as the angle between the jets in the plane 
perpendicular to the beam axis. A maximum value of 2.9 rad is important to 
reduce the remaining background. 

The result of this study is 4.3% detection efficiency and 9 background events. A 
detection confidence level of 3cr (99.73%) is expected for a cross section of 23 fb. 
Expected signal and background are shown in Fig. 18. 


Channel 

XiCXi c 

qq 

ww 

eWu 

tt 

ZZ 

eeZ 

Total (in 1000) 

1 

125 

70 

50 

7 

6 

60 

After Preselection 

391 

1652 

2163 

3185 

1259 

182 

318 

Evis/yfs < 0.4 

332 

202 

285 

3032 

70 

4 

98 

Njet = 2 

293 

172 

182 

2892 

17 

3 

72 

No isolated e or n 

218 

152 

98 

2757 

5 

3 

9 

E^ b /E vis < 0.3 

185 

101 

70 

2049 

5 

2 

4 

Invariant mass of jets>120GeV 

52 

25 

12 

7 

1 

0 

0 

Acoplanarity < 2.9rad 

43 

0 

5 

3 

1 

0 

0 


Table 3: Final event selection cuts, expected signal efficiencies, and the 
number of expected background events. Bold face numbers 
indicate major background reductions. 


The final xt bXi b event selection is summarized in Table |j. Here the cuts are: 

• A hard lower energy cut reduces most of the eWu background. 

• Topologies with back-to-back jets are reduced by an upper cut on the event 
thrust (Fig. 19). 

• A lower cut on the number of hadronic clusters reduces efficiently low- 
multiplicity background final states (Fig. 20). 

• Jets are clustered using the JADE algorithm. The y-cut value is optimized to 
obtain four jets for the signal. 

• Events with an isolated electron or muon are rejected. 

• An upper cut on the visible energy reduces qq, W + W~, and tt background. 

• Finally, the remaining tt background events are reduced by requiring less than 
30% perpendicular energy imbalance. 
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Concerning the number of 6-quarks per event, the decay XibXib —> W + XibW~Xib 
leads to the same final states as expected for tt background. Therefore, the tagging 
of 6-quarks has not proved to be efficient to reduce this background. 

The result of this study is 4.5% detection efficiency and 8 background events. A 
detection confidence level of 3a (99.73%) is expected for a cross section of 19 fb. 
Expected signal and background are shown in Fig. 21. 


Channel 

XibXib 

qq 

WW 

eSNv 

tt 

ZZ 

eeZ 

Total (in 1000) 

1 

125 

70 

50 

7 

6 

60 

After Preselection 

695 

1652 

2163 

3185 

1259 

182 

318 

Evis/Vs > 0.35 

610 

1494 

2011 

337 

1234 

178 

239 

Thrust <0.85 

536 

326 

420 

24 

1141 

69 

137 

Ncluster > 60 

399 

195 

134 

0 

769 

41 

3 

Njet = 4 

211 

53 

72 

0 

432 

22 

0 

No isolated e or /i 

99 

41 

49 

0 

105 

16 

0 

E vis /y/s < 0.55 

57 

3 

8 

0 

23 

0 

0 

Er b /E vis <0.3 

45 

1 

3 

0 

4 

0 

0 


Table 4: Final event selection cuts, expected signal efficiencies, and the 
number of expected background events. Bold face numbers 
indicate major background reductions. 


At a future e + e - collider with = 500 GeV, a large discovery potential for scalar 
top quarks is already expected within one year of data-taking (£ = 10 fb -1 ). Detector 
performances known from LEP detectors result in good background reduction. Full 
hermeticity of the detector is essential. 

The confidence levels for discovering a signal are shown in Figs. 22 a and 22 b for 
the xfbXib and Xi c Xib channels, respectively. Here, the confidence levels are given 
in a = Agxpected/ y/ ^background• The sensitivity is sufficient to discover a 200 GeV stop 


independently of the values of the mixing angle with 3a in both Xi c and X\b decay 
modes for the investigated neutralino and chargino mass combinations. A complete 
set of mass combinations remains to be studied. 

At a later stage, the total luminosity could reach 50 fb -1 , and the resulting discov¬ 
ery region for a 3 a effect is shown in Fig. 23. An increase of the center-of-mass en¬ 
ergy would extend the discovery region further as shown in Fig. 24 for yfs = 800 GeV 
C = 200 fb - . Based on experience made in the LEP2 searches |35j the efficiency 
for the simulated mass combination can be extended to a larger mass region for 
— m -± > 20 GeV. 
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6 Determination of Soft—Breaking Parameters 
A Case Study 


In this section we want to estimate the experimental accuracies for the stop and 
sbottom masses and mixing angles which can be expected from the Monte Carlo 
simulation described in the preceding sections. Without beam polarization a possible 
way to determine and cos Q~ t is using the yfs and cos 9f dependence of the unpo¬ 
larized e + e“ —> t\t\ total cross section (see Figs, la and 6 ). Let us take as reference 
point rn^ = 180 GeV, cos 6$ = 0.57, and y/s = 400 GeV and fs = 500 GeV as the 
two reference energies. Note that at | cos 6f = 0.57| the e + e _ —> t\t\ cross section 
has its minimum. The cross sections at this point for these two energy values are 
a = 18.2 ± 4.1 fb at fs = 400 GeV and a = 47.4 ± 5.5 fb at fs = 500 GeV where 
the experimental errors follow from the Monte Carlo simulation. Figure 25 shows 
the corresponding error bands in the rri^ — cos 6- t plane. As can be seen, hardly an 
information can be obtained on the mixing angle. 

The polarization of the e~ beam offers the possibility of measuring the sfermion 
masses and especially the mixing angles with much higher accuracy. The cross sections 
of e + e _ —> t\t\ for 90% left- and right-polarized e~ beam at the reference point = 
180 GeV, |cos0f| = 0.57 for fs = 500 GeV are a l = 48.6 ±6.0 fb, or = 46.1 ±4.9 fb, 
where the experimental errors are given by Acr/cr = N signal / N signal ± N background 
with the number of signal and background events determined as described in the 
previous section. Figure 26 shows the correponding error bands and the error ellipse 
in the — cos Of plane. The experimental accuracies obtained in this way for the 
mass of the lighter stop and the stop mixing angle are 

rn~ tl = 180 ±7 GeV, (27) 

cos 9 1 = 0.57 ±0.06. (28) 


We treat the sbottom system in an analogous way. Assuming that tan/3 is not 
too large we can neglect left-right mixing in the sbottom sector. In the “Minimal 
Super gravity-inspired Model” [|36| one expects m b < , thus bi = b L and b 2 = b r, 


i.e. cos 9~ b = 1. As reference point of the sbottom system we take = 200 GeV, 
m~ b2 = 220 GeV. The cross sections for e + e _ —> b\bi with 90% left-polarized e~ beams 


and for e + e —> b 2 b 2 with 90% right-polarized e beams then are crf,(e + e —> bibi) = 
61.1 ± 6.4 fb, <jR(e + e~ —> b- 2 b 2 ) = 6.0 ± 2.6 fb, where the errors are again determined 
by our Monte Carlo procedure. The errors for the sbottom masses follow as: 


m- bi = 200 ± 4 GeV, (29) 

m b2 = 220 ± 10 GeV. (30) 

With these values for rn ^, cos 6% m bi , and m b2 we can use (]T6|) and obtain the mass 
of the heavier stop t 2 if tan /3 is known from other experiments. Taking, for instance, 
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tan (5 = 2 leads to 


rn~ t2 = 289 ± 15 GeV. (31) 

Confirming this value by producing ^2 at higher energies would be an independent 
test of the MSSM. 

Assuming that also /i is known from an other experiment we are now able to 
calculate the underlying soft SUSY breaking parameters Mq, Mp, M b , A t and A b for 
the squarks of the third family according to eqs. (8) to (15). Taking /i = —200 GeV, 
tan (3 = 2, and m t = 175 GeV we obtain the following values: 


Mq 

= 195 ±4 GeV, 

(32) 

Mo 

= 138 ±26 GeV, 

(33) 

M b 

= 219 ±10 GeV, 

(34) 

A* 

= -236 ± 38 GeV if cos 0 t - > 0, 

(35) 

At 

= 36 ± 38 GeV if cos6^<0. 

(36) 


These results have to be compared with those of [Q and ||, where the stop mass 
was determined by a kinematical reconstruction of the bxt decay. In || a somewhat 
higher accuracy was obtained for unmixed squarks because one parameter less is 
involved in this case. 


7 Summary 

In this article we have discussed the production of stop, sbottom, r-sneutrino and 
stau pairs in e + e _ annihilation in the energy range = 500 GeV to 2 TeV. We 
have presented numerical predictions within the Minimal Supersymmetric Standard 
Model for the production cross sections and the decay rates and analyzed their SUSY 
parameter dependence. If tan/?> 10, not only the top Yukawa terms, but also the 
bottom and tau Yukawa terms have important effects. The production cross sections 
as well as the decay rates of stops, sbottoms and staus depend in a characteristic way 
on the mixing angles. 

A Monte Carlo study of e + e - —> t\ i\ at -y/s = 500 GeV with the decays t\ —> 
cxi and t\ —> bxt has been performed for m bl = 180 GeV, m^o = 100 GeV, and 
m bl = 180 GeV, m-± = 150 GeV, m^o = 60 GeV, respectively. A suitable set of 
kinematical cuts has been applied to reduce the known background reactions. In 
addition, detection regions have been given. 

We have also estimated the experimental accuracies for the masses of stops and 
sbottoms and the stop mixing angle from measurements of the polarized cross sections. 
Furthermore, we have made an estimate of the accuracies which can be obtained for 
the soft-breaking SUSY parameters. 

In summary, an e + e _ collider — especially with a polarized e~ beam — is an ideal 
machine for detecting and studying the scalar partners of the third generation quarks 
and leptons. 
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Figure Captions 


Fig. la: Contour lines for the total cross section of e + e —> t\ t\ in fb at a/s = 
500 GeV as a function of and cos Of . ( mi, = 300 GeV and rrig = 300 GeV.) 

Fig. lb: Total cross section of e + e —> t\ 1 1 in fb at y/s = 500 GeV as a function of 

cos 9f, for unpolarized (U) as well as left (L) and right (R) polarized e~ beams and 

rrif^ = 180 GeV. (m t ~ 2 = 300 GeV and rrig = 300 GeV.) 

Fig. 2a: Contour lines for the total cross section of e + e - —> £ 2^2 in fb at a/ s = 2 TeV 
as a function of mi 2 and cos 9f. (rn^ = 300 GeV and rrig = 700 GeV.) 

Fig. 2b: Total cross section of e + e —> ^£2 in fb at y/s = 2 TeV as a function of 

cos 9{, for unpolarized (U) as well as left (L) and right (R) polarized beams and 

rrif 2 = 700 GeV. (m bl = 300 GeV and rrig = 700 GeV.) 

Fig. 3a: Contour lines for the total cross section of e + e —> b\ b\ in fb at y/s = 
500 GeV as a function of m bi and cos 9- b . (mg = 300 GeV and rrig = 300 GeV.) 

Fig. 3b: Total cross section of e + e —> 61 61 in fb at y/s = 500 GeV as a function of 
cos $ 1 , for unpolarized (U) as well as left (L) and right (R) polarized beams and 
m~ bi = 180 GeV. (m b2 = 300 GeV and rrig = 300 GeV.) 

Fig. 4a: Contour lines for the total cross section of e + e~ —> f\ t 1 in fb at y/s = 
500 GeV as a function of rri^ and cos Of. 

Fig. 4b: Total cross section of e + e —> T\ t 1 in fb at y/s = 500 GeV as a function of 
cos 6 h, for unpolarized (U) as well as left (L) and right (R) polarized e~ beams and 
rrif x = 180 GeV. 

Fig. 5 : Total cross section of e + e —> u T i) r at y/s = 500 GeV as a function of 
for unpolarized (U) as well as left (L) and right (R) polarized e~ beams. 

Fig. 6: Total cross section of e + e —> i\ / at = 180 GeV and cos Q- t = 0.7 as a 
function of y/s. (rrif 2 = 300 GeV, rrig = 300 GeV.) 

Fig. 7: Gluon, gluino-top, initial state, and total radiative corrections relative to the 
tree level cross section of e + e~ —> G / at m bl = 180 GeV and cos% = 0.7 as a func¬ 
tion of y/s. (mf 2 = 300 GeV, rrig = 300 GeV.) 

Fig. 8a: Kinematically allowed parameter domains in the M — [i plane for rn^ = 
180 GeV and tan (3 = 2 for the decays: a) t\ —> c \?, b) t\ —> bW + Xi, c) t\ —> c 
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d) t\ —> bxt- e) ti —■> bxt- The grey area is covered by LEP2 for y/s = 190 GeV. 

Fig. 8b: Kinematically allowed parameter domains in the M — /j plane for mg = 

180 GeV and tan/? = 30 for the decays: a) b\ —> 6y5, b) 6 i 6y5, c ) 6 i ~^ 6 yg. The 
grey area is covered by LEP2 for y/s = 190 GeV. 

Fig. 9a: Branching ratios for the t\ decays as a function of M for m^ = 180 GeV, 
cos Of = 0.7, tan/3 = 2, and // = —500 GeV. The curves correspond to the following 
transitions: o —> cy^, ■ ti —► 6x1, —> 61F + y5- The grey area is covered by 

LEP2 for y/s = 190 GeV. 

Fig. 9b: Branching ratios for the b\ decays as a function of M for mg = 180 GeV, 
cos9- b = 0.7, tan/3 = 30, and /j = —500 GeV. The curves correspond to the following 
transitions: o b\ —> 6x5, D 6 i —> bx 2 . The grey area is covered by LEP2 for y/s = 190 
GeV. 

Fig. 10a: Branching ratios for the 6 i decays as a function of M for mg = 180 GeV, 
cos#g = 0.7, tan/3 = 30, and /i = —130 GeV. The curves correspond to the following 
transitions: o 6 1 —» Ay^, □ b± —> 6 ylj, A 6 i —> 6 x-j- The grey area is covered by LEP2 
for y/s = 190 GeV. 

Fig. 10b: Branching ratios for the f\ decays as a function of M for m^ = 180 GeV, 
cos 6f = 0.7, tan/3 = 30, and /i = —130 GeV. The curves correspond to the following 
transitions: o f\ —> tx?, □ Ti —> t x 2 , A fi —> r y°, O f\ —> ry^, ■ A —> ry x7- The 
grey area is covered by LEP2 for = 190 GeV. 

Fig. 10c: Branching ratios for the v y decays as a function of M for mp T = 180 GeV, 
tan/3 = 30, and /j = —130 GeV. The curves correspond to the following transitions: 
oj) T ^ Ur y°, □ P T —> /y r y°, y°, A 7y -> z/ T x°, ■ A -> t yl- The grey area 

is covered by LEP2 for y/s = 190 GeV. 

Fig. 11: Branching ratio for the decays of the heavier stop as a function of M for 
tan/3 = 2 and ji — —l TeV. 

(a) shows the decays into fermions: — t 2 —> t g, o t 2 —»• t y°, □ t 2 —>► t x 2 , ■ t 2 —> 6 yl. 

(b) shows the decays into bosons: ot 2 —> Zfi, □ t 2 —> 6 ° 61 , A t 2 —> H° t \, 
Ot 2 — 5 ► A°ti,u t 2 —> TV + 6 i. The grey area is covered by LEP2 for ^/s = 190 
GeV. 

Fig. 12: Branching ratio for the decays of the heavier sbottom as a function of M 
for tan/3 = 30 and fi = — 1 TeV. 

(a) shows the decays into fermions: — b 2 —■> 6 g, o b 2 —>■ 6 y 5 , □ b 2 —>► 6 y^, ■ b 2 —>► t yf. 

(b) shows the decays into bosons: o / 2 —> Z 6 lr □ b 2 —> 6 ° 61 , A b 2 —>► 77° 61 , 
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O & 2 — 1> A 0 b\, ■ b 2 —>► W ti, • b 2 —■► H ti. The grey area is covered by LEP2 
for y/s = 190 GeV. 

Fig. 13: Branching ratio for the decays of the heavier stau as a function of M for 
tan/3 = 30 and n = — 1 TeV. 

(a) shows the decays into fermions: o T 2 —> r Xi, □ fg —>► ■ fg —>► ^yyA 

(b) shows the decays into bosons: o —> Zfi, □ T 2 —>► h° f\ , A T 2 —>► H°f\, 
O T 2 — > A 0 fi. The grey area is covered by LEP2 for a/s = 190 GeV. 

Fig. 14: Branching ratio for the decays of the tau sneutrino as a function of M for 
tan/3 = 30 and n — — 1 TeV. 

(a) shows the decays into fermions: o z> T —>► v T Xi, uv T ^u T ■ f'r —> 

(b) shows the decays into bosons: u u T —> W + fi, • z> T —> H + f\. The grey area is 
covered by LEP2 for a/s = 190 GeV. 

Fig. 15: Background reactions and their cross sections for y/s = 500 GeV. 

Fig. 16: E vis /y/s < 0.4 for qq, WW, eWiq tt, ZZ, eeZ. 

Fig. 17: m inv > 120 GeV for x?cy?c, qq, WW, eWzq tt, ZZ, eeZ. 

Fig. 18: Sensitivity for an e + e“ —> t\t\ —> x^cy^c signal. Open histograms show the 
simulated signal, solid and hatched histograms show the remaining background after 
all selection cuts are applied. 

Fig. 19: Thrust < 0.85 for xt^Xib, qq, WW, eW u, tt, ZZ, eeZ. 

Fig. 20: Number of cluster > 60 for x^&Xi b , qq, WW, eW u, tt, ZZ, eeZ. 

Fig. 21: Sensitivity for an e + e“ —> t\ti —» XibXib signal. Open histograms show the 
simulated signal, solid and hatched histograms show the remaining background after 
all selection cuts are applied. 

Fig. 22: Detection confidence levels, (a) xtbXi b channel, (b) x'j’cy'^c channel. 

Fig. 23: Discovery region for yfs = 500 GeV and C = 50 fb -1 . In the shaded area a 
3 a effect is expected. 

Fig. 24: Discovery region for y/s = 800 GeV and £ = 200 fb -1 . In the shaded area 
a 3 a effect is expected. 

Fig. 25: Error bands for the total tree-level cross section of e + e —> tit\ in fb 
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at \fs = 400 GeV and yfs = 500 GeV as a function of rn^ and cos 6 j. The dot 
corresponds to rn^ = 180 GeV and cos Of = 0.57. The error bands are defined by 
( 0 - 400 , A 0 - 4 OO) = (18.2,4.1) fb and (cr 500 , Ao- 500 ) = (47.4,5.5) fb. 

Fig. 26: Error bands (dashed) and the corresponding error ellipse as a function 
of mfj and cos 9f for the total tree-level cross sections of e + e _ —> t\ti in fb at 
a /s = 500 GeV with 90% left- and right-polarized electron beam. The dot cor¬ 
responds to rn^ = 180 GeV and cos 61 = 0.57. The error bands are defined by 
(ffi, A(7 l) = (48.6,6.0) fb and (< 7 r,A< 7 r) = (46.1,4.9) fb. 
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